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Abstract 

We solve the non-relativistic Coulomb Shrodinger equation in d = 
2+1 via sine collocation. We get excellent convergence using a gener- 
alized sine basis set in position space. Since convergence in position 
space could not be obtained with more common numerical techniques, 
this result helps to corroborate the conjecture that the use of a local- 
ized basis set within the context of light cone quantization can yield 
much better convergence. All of the computations presented here 
were performed on an IBM-compatible PC with an Intel 486DX2-66 
microchip. 

Recently, light cone quantization (LCQ) of quantum field theory has at- 
tracted considerable attention as a possible alternative method for solving 
non-perturbative problems in quantum field theory In this scheme, one 

*preprint UTAH-IDR-CP-05 
T Email: koures@dirac.chem.utah.edu 



1 



obtains a rational, closed-form, and relativistically covariant Hamiltonian 
which allows one to avoid many of the severe mathematical difficulties which 
have plagued traditional equal time quantization techniques |]]. The discrete 
version of the light cone Hamiltonian allows one to model gauge theory on a 
computer as an eigenvalue and matrix diagonalization problem over a discrete 
and covariantly regularized Fock space0. 

A numerical implementation of the Discretized Light Cone Quantization 
(DLCQ) method to study positronium has produced promising results, but it 
is evident that the use of a plane wave basis will be a major bottleneck when 
one applies this method to more complicated models 0. It has been suggested 
that the use of a localized basis set would yield much better convergence for 
bound state computations 0, ||. The method introduced in Ref. [5] has now 
been adapted to the orthogonal localized basis set of Sine functions. 

Sine methods have increasingly been recognized as very powerful tools 
for attacking problems within applied physics and engineering ||. Until now, 
however, it was not clear how to apply Sine methods within the context 
of the operator formalism of quantum field theory. This general formalism 
will be reported in detail elsewhere; the focus of the present paper is to 
demonstrate the power of the sine collocation method by solving the radial 
Coulomb equation in d = 2+1. 

The motives for studying quantum electrodynamics (QED) in 2+1 di- 
mensions are numerous 0. The lower dimensions allow a smaller number 
of degrees of freedom but the model still possesses independent photon de- 
grees of freedom, unlike the (l+l)-dimensional model. The model is super- 
renormalizable and, when formulated with four-component spinors, it ex- 
hibits confinement 0, §. 



The non-relativistic Coulomb Schrodinger equation is derived from the 
LCQ formalism as follows @. One first derives the discretized light-cone 
Hamiltonian for (2+1) -dimensional QED with four- component spinors. A 
Tamm-Dancoff integral equation is then obtained for the "positronium" bound 
states. Taking the weak-coupling limit, one gets a non-relativistic integral 
equation which is the momentum-space Coulomb Schrodinger equation; the 
infrared divergences cancel between the self-mass and one-photon exchange 
diagrams [[?]]. The position-space result is then obtained by a Fourier trans- 
form, 

' 1 Q 2 

Vr +7r(l + lnmr ) $(r) = E$f(r) , (1) 

m Z7T 

where 7 is the Euler-Mascheroni constant, m is the mass, and g is the coupling 
constant. It is convenient to rewrite eq. ([!]) in terms of the dimensionless 
variables x = 



^f- r and A' 



+ In J^l. We get 
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or 



(- \J 2 X + lnx + 7 + In 2) V(x) = A'*(cc) 



with 



A = A' 



(2) 
(3) 
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The separation of variables, *ff(x) = R(x)Q(9), leads to 0(9) = exp(±iW), 
where I is the angular momentum quantum number. We are left with a radial 
differential equation for R which, after using the substitution H, 



R( x ) = x~ 1/2 f{x) 



transforms into a differential equation for /: 

"4/ 2 - 1 



f(x) + 



Ax 2 



+ lnx 



f(x) = Xf(x) 



(5) 



(6) 



This equation for / represents a singular Sturm-Liouville system and it can 
be solved by sine collocation||. 

To this end, we begin with the definition of the Sine function. If a > 
and m is an integer, the Sine function, S(m,a)(x), is defined by 

sin - (x — ma) 

S(m,a)(x) = — ±j r— L . (7) 

- [x — ma) 

The theory of sine series (cardinal functions) on the entire real line has been 
thoroughly developed. For a class of functions known as the Payley-Weiner 
class, the sine interpolation and quadrature formulas are exact |J. However, 
a more practical application of sine approximation of functions which are 
in a much less restrictive class has also been developed and the absolute 
errors have been derived via contour integration |J. In a nutshell, the Paley- 
Wiener functions are entire but a more practical class of functions should 
have specific growth restrictions on the real line and should be analytic only 
on an infinite strip centered about the real line, 

D s = {zeC : z = x + iy , | y | < d} . (8) 

C denotes the set of complex numbers. It turns out that the absolute er- 
ror of sine interpolation and quadrature on such functions is exponentially 
damped 0. 

For the problem of interest to us we need to use sine methods on a function 
/ whose domain is (0, oo). The more general class just described, however, 
has a domain which includes the whole real line. This conflict is elegantly 
circumvented through the use of conformal maps |J . For example, let be 
a one-to-one conformal map from some domain D to domain Ds and let ip 
denote the inverse map which is also conformal. If / is analytic in D then 



/ o ip is analytic in D$- So, if a numerical process has been developed in a 
domain containing the whole real line, 1Z, then this process can be carried 
over to a new domain containing only a proper subset of the real line. 
We will shortly see that the domain of interest to us is 

D = ju?eC : | axg(sinh(w)) | < d < . (9) 

This domain is conformally mapped onto the infinite strip D$ by the function |i 

z = 4>{w) = \n(sinh(w)) . (10) 

If we let 



w = i/j( z ) = <jr\z) = \n(e z + Vl + e 2z ) , (11) 

then 

r = = (o,oo) , (12) 

as desired. For a > and m an integer we define the sine points 

x m = ^(mo) = In (e ma + V 1 + e 2ma ) . (13) 
One may further verify that 



^'(x m ) = VI + e" 2ma 



x m ) = -e~ 2ma , (14) 

where the primes denote differentiation with respect to x. 

Given the above definitions and maps, suppose we have a function / and 
that there are positive constants a, (3, and C such that 

. x a , xe fo,ln( 1 + y/2)) , N 

e p:r , xe ln(l + V2),oo) 



If we choose 



N = ceil{^M) 

P 



V aM ~ ln(2) 
where ceil{x) rounds x to the nearest integer > x, then 

A' 



2nd 2nd . . 

< 77^ » ( 16 ) 



f(*)- E /(x m )S(m,a)(0(a;)) |= O ( exp - V ndaM ) , (17) 



m=-M 

and 



— Yl f(x m )3{m,a)(<l>(x)) \= 0(exp-VndaM . (18) 



dx k dx k m= _ M 

As promised, we see that the absolute errors are exponentially damped. 

The two equations ( fL7|) and (p~8|) are all we need to solve the differential 
equation (|B|). In addition, we would like to have a normalized final answer 
and so we also introduce the result for sine quadrature ||. Suppose we have 
a function F and that there are positive constants a, (3, and C such that 



F(x) I < C 



Thenl 



I [°°F(x)dx-a Y ^4 1= Cfexp-v / 27rdaM N ) , (20) 

with M, N, and a being defined the same as before. 

Now, to solve (||) we note that as x — > the physically acceptable solution 
takes the form[§] 

/ oc x l+1 ' 2 . (21) 



(o,ln(l + v^) 
ln(l + V2),oo 



(19) 



From ( jl5|) we may thus take a = I + 1/2. On the other hand, as x — ■> oo 
the physically acceptable solution takes the form 



/ oc e 



-x In a; 



(22) 



So, from ( |T5|) we may take /3 < 1. Next, we use ( |17|) and CflSf) to approximate 
/ and /", respectively: 



/(*) 

fix) « v 



^ S(m,a)((f>(x))f(x m ) 
d 2 



m=—M 
N i 



(</>'(*)) 2 



(23) 

S(m,a)(0(x)))/(x m ).(24) 



m=-M 

To evaluate these expansions at a general sine point x n note that 



SV,a)(0(z))k = <L 

— S{m,a){(f){x)) \ Xn = -6$ m 
deb a 



d 2 



S(m,a)((f)(x)) \ Xn 



L 5 m 

a 2 n ' m 



where 5$' m is the Kronecker delta function and 
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"n,m 

si 2 l = 
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Using these results along with (I|) we may now approximate 
collocation, 

N 



[-^n.m + ^n,m -^n,m f( x m) — ^f( x 



(25) 



(26) 
(27) 



via smc 



(28) 



where 



r(o) 

n,m 
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In the tradition of LCQ, we have formulated the problem as an eigenvalue 
and matrix diagonalization problem. The components of the eigenvector 
f(x m ) can be substituted into (p3~l) to compute the eigenfunction at arbitrary 
x. Finally, we use (§) to acquire R(x) and we normalize this function by 



computing the norm, J R 2 (x) dx , using the sine quadrature formula 

In Table 1 we list the first five eigenvalues for I ranging from to 4. The 
diagonalization was performed with MATL AB |9|] . We used d = 7r/4, (3 = 
1/2 to 1, and values for M up to 500. The convergence was excellent for all / 

Table 1: First five eigenvalues for I ranging from to 4. 





/ = 


I = 1 


I = 2 


I = 3 


I 4 


A 


0.52643626 


1.3861862 


1.8443720 


2.1578468 


2.3962798 


Ai 


1.6619365 


2.0094748 


2.2758614 


2.4881158 


2.6638815 


A 2 


2.1870578 


2.3943387 


2.5800522 


2.7390550 


2.8772701 


A 3 


2.5153639 


2.6726676 


2.8144703 


2.9409664 


3.0543788 


A 4 


2.7677810 


2.8906069 


3.0049630 


3.1096821 


3.3373990 



values but we needed considerably higher M values for / = than for I ^ 0. 
In Ref. [7] the first five / = eigenvalues were computed in momentum space: 

Aq = 1.7969 , A; = 2.9316 , A^ = 3.4475 

A' 3 = 3.7858 , A^ = 4.0380 , (30) 



where A' is related to A in Table 1 by (f|) . These momentum space results do 
not agree with the position space computations f|] and additional repeated 
attempts to get convergence in position space have failed^. It was concluded 
in Ref. [7] that "the previous position space calculation || was inaccurate, due 
to the slow, logarithmic behavior of the potential . . . The momentum space 
calculation is much more rapidly convergent" . 

On the other hand, if we use (^) to convert the first column of Table 1 to 
the parameterization used in Ref. [7] we see that 

X' Q = 1.7967991 , X[ = 2.9322993 , A' 2 = 3.4574206 

A3 = 3.7857268 , A^ = 4.0381439 ; (31) 

These values are in very good agreement with (|3"0"|), they converge faster, 
and they are more accurate. Note further that the eigenvalues for / 7^ are 
actually in good agreement with the position space results M. (The momen- 
tum space calculations in [7] were not carried out for Z 7^ 0.) This leads us 
to conclude that the slow convergence of the previous position space results 
was not due to the slow logarithmic behavior of the potential; rather, it was 
due to a minor instability caused by the sign-flip of the "centrifugal" part of 
the potential when 1 = 0. 

As promised, it is now very easy to use (^) along with @ to compute 
any eigenfunction; these eigenfunctions are easily normalized through the use 
of the sine quadrature result (|20|) . As a representative sample, in Figures (1) 
to (3) we display the first three normalized eigenfunctions for I = 0; note 
that the number of bumps increases incrementally as we go from n = to 
n = 2. Figures (4) to (6) show the first three normalized eigenfunctions for 
1 = 4; note that the "centrifugal" barrier causes the bumps to move away 



from the center, relative to the I = plots. 

We have thus demonstrated that sine methods provide a very powerful 
tool for solving the radial Coulomb Schrodinger equation in d — 2 + 1. These 
methods are very accurate and converge very fast. This point is strongly 
demonstrated by the fact that convergence for the / = values could not be 
attained in position space by other well-known numerical techniques. The 
sine collocation method used here may be straightforwardly extended to at- 
tack more complicated problems within quantum field theory; this formalism 
will be reported in detail in the future. 
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Figure 1: Plot of normalized function Ro,o(x) for / = and A = 0.5264. 



Figure 2: Plot of normalized function i2 lj0 (x) for I = and Ai = 1.6619. 



Figure 3: Plot of normalized function R,2,o(x) for / = and A 2 = 2.1871. 



Figure 4: Plot of normalized function Rq^x) for / = 4 and A = 2.3963. 



Figure 5: Plot of normalized function R\^{x) for / = 4 and Ai = 2.6639. 



Figure 6: Plot of normalized function R2,±(x) for / = 4 and A 2 = 2.8773. 
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